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Acoustic  backscatter  returns  measured  with  a  vertical  source  and  horizontal  receiving  array  are 
beamformed  and  charted  to  respective  scattering  sites  on  the  western  flank  of  the  Mid-Atlantic 
Ridge  (MAR).  Well-defined  patterns  of  high-level  backscatter  resembling  bottom  morphology 
occur  within  the  direct  path  and  at  convergence  zone  ranges.  Measured  backscatter  is  compared 
with  modeled  two-way  transmission  loss  using  the  wide-angle  parabolic  equation  (PE)  and 
high-resolution  supporting  bathymetry.  The  level  of  backscatter  is  found  to  be  strongly 
dependent  upon  two-way  transmission  loss  (TL).  Specifically,  TL  selects  prominent 
bathymetric  features  from  which  high  backscatter  is  returned,  with  backscatter  maxima 
corresponding  to  TL  minima.  Scattering  strength  estimated  for  the  wide  area  surveyed  shows 
low  variance  and  relatively  minimal  spatial  variation.  Apparently,  distinct  scattering 
contributions  are  so  well  mixed  in  long-range  backscatter,  that  dominant  variations  are  generally 
due  to  propagation.  These  results  indicate  that  it  is  possible  to  forecast  the  range  and  azimuth 
of  prominent  MAR  backscatter  returns  simply  by  propagation  modeling  if  high-resolution 
bathymetry  is  available.  They  also  have  led  to  the  development  of  a  new  technique  for  ambiguity 
removal  in  reverberation  data  measured  with  a  horizontal  line  array.  The  method  takes 
advantage  of  environmental  asymmetry  provided  by  TL  in  highly  range-dependent 
environments  and  does  not  require  multiple  measurements  as  do  other  existing  methods. 
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INTRODUCTION 

In  July  and  August  of  1991  acoustic  reverberation 
data  were  acquired  above  the  western  flank  of  the  Mid- 
Atlantic  Ridge  (MAR).  A  single  research  vessel  towed  a 
vertical  source  and  horizontal  receiving  array  to  obtain 
quasimonostatic  measurements  of  ocean-bottom  reverber¬ 
ation  at  low  frequency.  The  three-week  experiment  covered 
a  large  area  and  employed  a  broad  range  of  waveforms  as 
well  as  impulsive  charges.  Data  were  recorded  long  enough 
to  address  a  variety  of  direct-path  and  long-range  reverber¬ 
ation  issues.  Since  sea  states  were  generally  below  Beaufort 
Force  1,  interference  from  surface  reverberation  was  neg¬ 
ligible  and  data  quality  was  high.  Referred  to  as  the  Acous¬ 
tic  Reconnaissance  Cruise  (ARC),  the  experiment  is  part 
of  an  on  going  Special  Research  Program  (SRP)  on  bot¬ 
tom  reverberation  sponsored  by  the  Office  of  Naval 
Research.1,2  In  August  of  1992,  a  supporting  database  of 
high-resolution  bathymetry3  was  obtained  for  the  region 
previously  surveyed  by  the  ARC. 

Our  goal  is  to  understand  the  relationship  between 
backscatter  returning  from  multiple  convergence  zone 
(CZ)  ranges  and  complex  geomorphology  in  the  deep 
ocean.  Previous  investigators  have  found  a  correspondence 
between  long-range  backscatter  measured  via  towed  array 
and  gross  bathymetry.4"8  In  an  analysis  of  Tyrrhenian  sea 
data,  Nghiem-Phu  and  Tappert  have  shown  that  range- 
dependent  modeling  of  two-way  transmission  loss  (TL) 
with  the  parabolic  equation  (PE)  can  determine  the  gen¬ 
eral  characteristics  of  backscatter  without  detailed  infor¬ 
mation  about  individual  scattering  sites.9  However,  this  as 
well  as  other  past  studies  have  focused  upon  environments 
where  isolated  large-scale  scattering  features,  such  as  sea¬ 


mounts  or  continental  margins,  penetrate  or  curtail  sound 
channels  which  are  otherwise  generally  free  from  bottom 
interaction.  In  these  situations,  gross  bathymetric  surveys 
are  sufficient  for  a  first-order  modeling  of  TL. 

The  ARC  stands  apart  because  it  was  carried  out  in 
the  previously  uninvestigated  MAR  where  rugged  bathym¬ 
etry  projects  small-scale  features  into  the  deep  sound  chan¬ 
nel.  As  a  result,  bottom  limitation  of  deep-sound-channel 
propagation  is  extensive  and  highly  variable  for  sources  at 
depths  used  in  typical  sonar  systems.  Reverberation  in 
such  complex  environments  has  not  been  thoroughly  in¬ 
vestigated  in  the  past,  and  is  interesting  because  of  the 
numerous  returns  found  over  range  and  azimuth  which 
generally  serve  as  false  targets  in  active  sonar  measure¬ 
ments.  Only  with  the  recent  availability  of  high-resolution 
bathymetry,  collected  in  support  of  SRP  acoustic  data,  has 
it  become  possible  to  study  the  relationship  between  rever¬ 
beration  and  small-scale  geomorphology  in  the  MAR. 
Such  analysis  is  not  only  useful  for  sonar  design  in  highly 
reverberant  environments,  but  may  also  provide  a  practical 
method  for  remote  survey  of  complex  bathymetry.  In  this 
paper,  we  show  that  long-range  backscatter  is  highly  cor¬ 
related  with  small-scale  geomorphology  and  that  this  cor¬ 
relation  can  be  largely  explained  by  TL  modeling  with 
high-resolution  bathymetry. 

To  describe  our  method  of  charting  backscatter  we 
find  it  convenient  to  introduce  the  following  definitions. 
We  define  an  “observation”  as  a  single  active  insonification 
and  subsequent  towed-array  measurement  of  reverberation 
where  the  array  position  and  orientation  can  be  considered 
constant.  We  define  a  “reverberation  map”  as  the  magni¬ 
tude  of  acoustic  returns  from  a  single  observation  charted 
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to  the  horizontal  location  of  respective  scattering  sites.  To 
generate  reverberation  maps  from  ARC  data,  we  first  esti¬ 
mate  the  range  and  azimuth  of  scattering  sites  by  two-way 
travel  time  analysis  and  beamforming.  Since  beamformed 
line-array  data  have  an  inherent  right-left  ambiguity  about 
the  array’s  axis  it  is  especially  challenging  to  identify  true 
from  virtual  scattering  sites  in  the  MAR  where  ambiguous 
returns  are  broadly  distributed  in  range  and  azimuth. 

The  ARC  tow-ship  tracks  were  designed  to  provide  a 
general  survey  of  backscatter  in  the  vicinity  of  the  western 
MAR.  They  are  not  optimal  for  the  existing  methods  of 
ambiguity  removal  which  rely  upon  multiple  adjacent  ob¬ 
servations  at  different  array  heading.10,11  Therefore,  we  in¬ 
troduce  a  method  based  upon  two-way  TL  modeling  with 
the  wide-angle  PE,12  and  use  this  in  our  subsequent  data 
analysis.  This  method  has  a  strong  physical  foundation 
which  takes  advantage  of  naturally  occurring  asymmetries 
in  the  MAR  environment  by  incorporating  the  high- 
resolution  bathymetry  database.  Unlike  other  ambiguity 
resolution  methods,  it  does  not  require  multiple  data  ob¬ 
servations.  Further,  it  can  be  applied  objectively  to  the 
entire  data  set.  This  is  in  contrast  to  ambiguity  resolution 
by  comparison  with  bathymetry  which  is  extremely  subjec¬ 
tive,  and  biased  toward  discrete  or  well-defined  returns. 

Once  ambiguity  is  resolved,  backscatter  levels  are  un¬ 
ambiguously  projected  onto  bathymetry  for  comparison. 
We  find  that  prominent  returns  occur  at  n-\-\  CZ  ranges 
and  are  well  correlated  with  ridges  which  have  extended 
scarps  facing  the  observation,  where  n= 0,1  for  the  present 
survey.  However,  resolution  is  too  poor  to  determine 
whether  this  correlation  is  solely  due  to  a  propagation  ef¬ 
fect  which  insonifies  the  backfacing  scarps  of  prominent 
ridges,  or  whether  it  is  also  due  to  the  orientation  and 
geomorphology  of  scatterers  within  the  ridges.  Finally,  we 
combine  modeled  TL  with  unambiguous  backscatter  level 
to  estimate  scattering  strength  over  the  wide-area  sur¬ 
veyed.  Our  technique  introduces  two  advances  over  more 
traditional  scattering  strength  estimation:  (a)  TL  is  com¬ 
puted  using  a  full-field  propagation  model  rather  than  a 
raytrace;  (b)  the  area  term  of  the  sonar  equation  is  incor¬ 
porated  by  a  spatial  convolution  rather  than  a  constant 
scale  factor. 


I.  A  HIGHLY  BOTTOM-LIMITED  ENVIRONMENT 
A.  Observation  geometry 

We  analyze  only  a  small  fraction  of  the  total  ARC 
reverberation  data  set  for  the  present  study,  focusing  on 
data  measured  central  to  the  bathymetry  database.  We 
classify  observations  by  track,  or  run  number  (R),  and 
emission  sequence,  or  epoch  number  (E).  The  observa¬ 
tions  for  the  present  study  are  then  labeled  R7E10  and 
R7E15  in  Fig.  1  where  track  information  is  superposed  on 
an  image  generated  from  the  high-resolution  bathymetry, 
sampled  at  200-m  intervals.  For  example,  track  7  runs 
from  north  to  south. 
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B.  Long-range  propagation  corridors 

Some  general  characteristics  of  the  acoustic  field  inci¬ 
dent  upon  the  ocean  bottom  can  be  readily  deduced  from 
source,  sound  speed,  and  bathymetric  information.  The 
source/receiver  deployment  geometry  is  indicated  in  Fig. 
2,  along  with  miscellaneous  operational  characteristics  of 
the  arrays.  Due  to  negligible  horizontal  variation,13  a  single 
representative  sound-speed  profile,  shown  in  Fig.  3,  can  be 
used  to  describe  the  water  column.  While  the  density, 
sound  speed,  and  attenuation  of  the  bottom  vary  in  both 
the  horizontal  and  vertical,  information  about  these  pa¬ 
rameters  is  limited.  However,  we  do  know  that  sediment  of 
varying  thickness  and  relatively  homogeneous  acoustic  im¬ 
pedance  covers  a  majority  of  the  survey  area.14  Therefore, 
we  assume  the  bottom  is  comprised  of  sediment  with  uni¬ 
form  density  1.5  g/cm3,  sound  speed  1700  m/s,  and  atten¬ 
uation  0.5  dB/A  for  all  range-dependent  acoustic  modeling 
done  in  this  paper. 

An  example  of  deep  sound  channel  propagation  char¬ 
acteristic  to  the  ARC  is  shown  in  Fig.  4  for  a  range-depth 
cross  section  of  the  ocean.  An  idealized  flat  bottom  at  5-km 
depth  is  used  to  demonstrate  the  CZ  structure  of  the  source 
beam  when  undisturbed  by  bottom  interaction,  which  is  a 
rare  situation  for  our  survey.  Propagation  modeling  is  done 
with  the  wide-angle  PE  for  the  environmental  parameters 
discussed  above.  The  source  array  is  accurately  modeled  by 
ten  omnidirectional  monopoles  at  depths  corresponding  to 
those  of  the  actual  projectors.  The  model  source  elements 
are  horizontally  steered  and  operate  at  270  Hz.  The  com¬ 
bined  TL  is  normalized  so  that  a  single  source  level  may  be 
applied  to  obtain  sound-pressure  level  at  any  point  in  the 
field. 

The  source  beam  contains  the  highest  amplitude  sound 
waves  by  orders  of  magnitude.  The  beam  has  a  complex 
structure  due  to  interference  caused  by  the  source’s  prox¬ 
imity  to  the  pressure  release  surface;  the  numerous  side- 
lobes  of  the  untapered  source  array;  and  the  refractive 
properties  of  the  water  column  which  lead  to  deep-sound- 
channel  propagation.  It  has  a  first  deep  vertex  approxi¬ 
mately  33  km  from  the  source,  which  defines  the  \  CZ 
range,  the  region  within  which  is  referred  to  as  the  direct- 
path  area.  Here  beam  vertex  occurs  between  roughly  3900- 
and  4500-m  depth.  For  practical  considerations  then,  the 
beam  will  propagate  undisturbed  in  the  water  column  for 
bathymetry  exceeding  4500  m,  will  be  partially  blocked  or 
reflected  for  bathymetry  between  4500  and  3900  m,  and 
will  be  totally  blocked  or  reflected  for  bathymetry  shal¬ 
lower  than  3900  m. 

To  consolidate  this  information,  we  employ  the  follow¬ 
ing  terminology.  “Conjugate  depth”  refers  to  the  depth 
below  the  minimum  sound  speed  where  the  sound  speed 
returns  to  the  value  at  the  source  center.  Conjugate  depth 
occurs  at  roughly  3900  m  for  the  ARC.  (By  Snell’s  law, 
horizontal  rays  emanating  from  the  source  will  vertex  at 
conjugate  depth.)  The  term  “excess  depth”  refers  to  the 
difference  between  bathymetric  depth  and  conjugate  depth. 
And,  “positive  excess  depth”  occurs  for  bathymetric 
depths  exceeding  the  conjugate  depth. 

To  give  an  overview  of  the  anticipated  blockages  in 
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FIG.  1.  Hydrosweep  bathymetry  from  SRP  geophysical  survey  of  1992  gridded  at  200-m  intervals.3  The  dotted  white  line  indicates  the  tracks,  or  straight 
line  runs,  of  the  research  vessel.  The  colored  dots  indicate  the  position  of  monostatic  backscatter  observations  relevant  to  this  investigation.  The 
observations  are  labeled  R7E10  and  R7E15  for  the  northern  and  southern  observations,  respectively.  The  irregular  border  within  the  figure  indicates  the 
limit  of  hydrosweep  coverage. 


multiple  CZ  propagation,  we  have  plotted  the  bathymetry 
where  positive  excess  depth  is  available,  in  Fig.  5,  and  note 
that  this  comprises  less  than  50%  of  the  area  surveyed. 
Two  concentric  circles  centered  about  observation  R7E15, 
closest  to  the  track  crossing,  are  superposed  on  the 


Research  vessel 


Tilt  <  2° 


vertical  Source  Array 

FIG.  2.  A  sketch  of  the  research  vessel  towing  the  vertical  source  and 
horizontal  receiving  arrays.  Array  depths  and  orientation  change  with 
ship’s  speed.  Conditions  shown  are  for  a  typical  tow  speed  of  3.5  kn.  The 
element  spacing  is  3.66  m  for  the  vertical  array  and  2.50  m  for  the  hor¬ 
izontal  array. 
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bathymetry.  These  indicate  the  \  and  CZ  radii  as  if  there 
was  no  bottom  interaction.  In  reality,  however,  bottom 
interaction  is  extensive.  The  \  CZ  circle’s  intersection  with 
positive  excess  depth  is  limited  to  relatively  narrow  azi¬ 
muths  to  the  southwest  and  southeast  of  the  observation. 
Only  along  these  restricted  azimuths  do  we  anticipate  mul¬ 
tiple  CZ  propagation  to  be  possible.  At  all  other  azimuths, 
bottom  interaction  within  the  direct-path  area  will  cause 
the  source  beam  to  be  attenuated  and  redistribution  in  ver¬ 
tical  angle,  resulting  in  much  weaker  long-range  insonifi- 
cation.  Therefore,  only  a  sparse  azimuthal  distribution  of 
scattering  sites  at  CZ  range  will  be  insonified  by  water¬ 
borne  paths.  This  is  significant  because  the  strongest  back¬ 
scatter  is  expected  to  return  from  sites  where  waterborne 
paths  are  available. 

To  demonstrate  the  environmental  symmetry  breaking 
potential  available  in  range-dependent  propagation  model¬ 
ing,  we  have  strategically  selected  radials  for  PE  runs.  The 
radials  extend  from  observation  position  R7E15  and  are 
comprised  of  ambiguous  pairs  (a, a')  and  (£,/?')  about  the 
array  axis.  Radials  a  and  0,  respectively,  intersect  the 
southwestern  and  southeastern  azimuths  where  positive  ex¬ 
cess  depth  is  available  at  \  CZ.  This  is  displayed  in  Fig.  5. 
Bathymetry  from  the  high-resolution  database  is  extracted 
along  each  radial  to  supply  range  dependence  in  each  run. 
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FIG.  3.  A  sound-speed  profile  typical  of  those  measured  via  XBT  during 
the  ARC.  The  cast  was  taken  approximately  18-km  south  of  the  central 
observation  R7E15  on  run  7.  Archival  data  used  below  the  2000-m  depth. 

To  be  consistent  with  ARC  steering  and  source  frequency 
for  the  data  analyzed  in  this  paper,  the  modeled  source 
elements  are  steered  to  a  depression  angle  6°  from  horizon¬ 
tal  and  operate  at  270  Hz.  Range-depth  cross  sections  of 
the  propagation  along  each  radial  are  shown  in  Fig.  6. 
First,  consider  propagation  along  radial  a  where  a  water¬ 
borne  path  is  available  to  a  bottom  site  at  roughly  \\  CZ. 
Here  TL  is  substantially  lower  than  at  the  ambiguous  \\ 
CZ  site  along  radial  a'.  This  is  because  along  a',  the  source 
beam  is  terminated  within  the  direct-path  area.  Therefore, 
backscatter  of  much  higher  amplitude  is  expected  from  the 
1^  CZ  bottom  site  along  a.  Propagation  runs  for  ft  vs  /?',  in 
Fig.  6,  show  a  similar  scenario.  That  is,  a  waterborne  in- 
sonification  of  the  bottom  at  1^  CZ  occurs  along  ft  but  not 
along  p\  [One  additional  observation  useful  to  understand- 
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FIG.  4.  Transmission  loss  (TL)  for  a  range-independent  cross  section  of 
the  ocean  using  the  sound-speed  profile  of  Fig.  3  and  the  ten-element 
vertical  source  array  shown  in  Fig.  2.  Convergence  zone  structure  typical 
of  the  region  is  evident  as  are  sidelobes  in  the  beampattem  due  to  the 
untapered  array,  which  is  steered  to  0°  depression  angle  for  this  example. 
Modeling  was  done  with  the  wide-angle  PE.12 


FIG.  5.  Positive  excess  depth  with  respect  to  the  conjugate  depth  of  the 
source-array  center.  The  white  dotted  lines  indicate  the  approximate  lo¬ 
cation  of  the  source-beam  vertex  at  j  and  lj  CZ  ranges  from  central 
observation  R7E15  which  is  indicated  by  a  dark  dot.  The  solid  white  lines 
indicate  radial  pairs  (a,a')  and  (ft/?')  symmetric  about  the  array  axis. 

ing  the  backscatter  process  should  be  made.  For  each  of 
the  (a, a')  and  (ft, ft')  runs  presented  in  Fig.  6,  backfacing 
scarps  are  generally  insonified  by  the  source  beam  while 
forward-facing  scarps  are  in  the  acoustic  shadow.  Stronger 
backscatter  is  therefore  expected  from  backfacing  scarps.] 

These  examples  illustrate  the  substantial  differences  in 
TL  that  can  exist  at  ambiguous  scattering  sites  due  solely 
to  asymmetry  in  bathymetry  along  ambiguous  radials.  In 
later  sections  we  extend  these  results  to  objectively  elimi¬ 
nate  right-left  ambiguity  from  backscatter  data  over  wide 
areas  by  comparing  TL  at  ambiguous  scattering  sites. 

II.  REVERBERATION  MAPS 
A.  Data  acquisition  and  reduction 

Although  a  wide  variety  of  waveforms  were  used  dur¬ 
ing  the  ARC,  from  narrow-band  continuous  wave  (cw) 
tonals  to  broadband  hyperbolic  frequency-modulated 
(hfm)  emissions,  we  limit  our  present  analysis  to  backscat¬ 
ter  from  cw  signals  of  r=2-s  pulse  duration.  This  is  prin¬ 
cipally  because  the  corresponding  range  resolution  A r 
—ct/ 2— 1.5  km,  for  a  nominal  mean  sound  speed  of  c=1.5 
km/s,  is  on  the  order  of  the  cross-range  resolution  for  the 
ranges  relevant  to  the  present  study.  (Cross-range  resolu¬ 
tion  is  simply  the  azimuthal  resolution  of  the  horizontal 
receiving  array  multiplied  by  the  range  of  the  scattering 
site.) 

Azimuthal  dependencies  of  the  backscatter  data  are 
determined  by  beamforming  time  series  data  received  on 
the  array  elements.  Data  were  beamformed,  filtered,  and 
demodulated  to  a  128-Hz  sample  rate  from  the  original 
rate  of  2048  Hz.  The  data  were  then  refiltered  to  within  a 
3-Hz  band  about  the  270-Hz  center  frequency  to  eliminate 
out-of-band  noise,  leading  to  the  beam-time  format  neces¬ 
sary  for  our  mapping  procedure. 
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FIG.  6.  Transmission  loss  (TL)  for  range-dependent  cross  sections  of  the  ocean  along  the  radials  depicted  in  Fig.  5.  In  each  case,  propagation  to  l|  CZ 
is  partially  or  fully  blocked  in  the  direct  path  for  one  member  of  each  symmetric  pair,  and  is  clear  for  the  other.  (The  source  is  modeled  with  the 
wide-angle  PE  as  an  untapered  ten-element  vertical  array,  shown  in  Fig.  2,  steered  to  a  6°  depression  angle.)  In  general,  backfacing  scarps  are  insonified 
by  the  source  beam  while  forward-facing  scraps  are  in  the  acoustic  shadow.  Stronger  backscatter  is  therefore  expected  from  backfacing  scarps. 
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FIG.  7.  A  reverberation  map  for  observation  R7E15.  Highest  level  returns  occur  within  the  direct  path  and  at  \\  CZ  range.  All  returns  are  mirrored 
about  the  array  axis,  roughly  183°  azimuth,  due  to  right-left  ambiguity.  Poor  resolution  due  to  endfire  occurs  within  roughly  ±30°  of  the  array  axis. 


scatter,  bathymetry,  and  TL.  Referring  to  Fig.  6,  backfac- 
ing  scarps  are  insonified  by  waterborne  paths  at  both  sites 
a  and  j8,  while  forward-facing  scarps  are  in  the  acoustic 
shadow.  The  scarp  at  site  a,  where  a  broad  backscatter 
return  is  charted,  corresponds  to  the  face  of  a  prominent 
abyssal  hill,  while  the  scarps  at  p,  where  narrower  returns 
are  charted,  belong  to  smaller  ridges  distributed  on  a  shal¬ 
lower  ledge. 


FIG.  8.  A  histogram  of  wide-area  backscatter  levels  for  observation 
R7E15  with  mean  level  73  dB  re:  1  /zPa,  median  72  dB  re:  1  /iPa,  and 
standard  deviation  8.5  dB.  Correcting  for  array  gain  and  finite  bandwidth, 
the  noise  floor  corresponds  to  typical  deep-ocean  ambient  noise  levels. 
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III.  TRANSMISSION  LOSS  MAPS 

For  comparison  with  reverberation,  we  are  interested 
in  the  two-way  TL,  which  is  proportional  to  reverberation 
level  R  in  the  sonar  equation 

R  =  -TLf-TLb+S+  10  log  A+W,  (3) 

where  S  is  defined  as  scattering  strength,  W  is  the  source 
strength,  A=rAr(3(0)  is  the  scattering  area,  r  is  the  range 
to  the  center  of  A ,  TL^  is  the  forward  transmission  loss 
from  the  vertical  source  array  to  bottom  sites,  and  TL*  is 
the  back  transmission  loss  from  bottom  sites  to  the  hori¬ 
zontal  receiving  array. 

To  map  TL y  and  TL*  over  the  wide  areas  surveyed  in 
each  backscatter  observation,  we  extract  bathymetry  along 
R  =  1 10-km  radials  at  evenly  spaced  azimuthal  increments 
A0=1°  about  each  observation.  Along  each  radial,  we 
compute  the  acoustic  field  in  a  two-dimensional  range- 
depth  cross  section  of  the  ocean  using  the  wide-angle  PE 
with  the  same  frequency  and  environmental  parameters 
used  in  Sec.  I,  Fig.  6.  The  TL  at  1  m  above  the  bottom 
(roughly  A/6)  is  extracted  at  evenly  spaced  ranges  along 
each  radial.  For  each  observation,  the  radials  are  combined 
to  produce  a  TL  map  over  the  entire  survey  area.  To  fill  the 
areas  between  radials,  the  TL  within  +/— A0/2  of  each 
radial  is  given  the  value  at  the  radial  without  interpolation. 
For  TL y  we  use  the  same  vertical  source  array,  steering, 
and  normalization  as  in  Fig.  6.  For  TL*  we  use  a  point 
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FIG.  9.  Forward  and  back  transmission  loss  maps  TLy  and  TL*  for 
observation  R7E15.  Modeling  is  done  for  radials  at  1°  increments.  A 
general  correspondence  can  be  found  between  TL  minima  in  these  figures 
and  backscatter  maxima  in  Fig.  7. 

source  at  the  depth  and  central  range  of  the  horizontal 
receiving  array  and  invoke  reciprocity,  with  the  under¬ 
standing  that  variations  in  TL*  along  the  receiving  array’s 
aperture  are  negligible.  Maps  of  TL f  and  TL*  are  shown  in 
Fig.  9.  While  there  is  a  general  correspondence  between  the 
two  maps,  sharper  features  are  found  in  TLy  due  to  the 
narrow  beamwidth  of  the  vertical  array.  This  is  especially 
so  in  the  direct-path  region  where  the  incident  field  of  the 
point  source  is  more  uniformly  distributed  on  the  ocean 
bottom. 

The  TLy+TL*  map  for  observation  R7E15  is  shown 
in  Fig.  10.  Direct-path  and  long-range  sites  insonified  by 
waterborne  paths  have  the  lowest  two-way  TL.  These  oc¬ 
cur  within  \  and  at  \\  CZ  ranges,  respectively,  the  latter 
being  contingent  upon  the  availability  of  a  long-range 
propagation  corridor.  Two-way  TL  to  sites  at  intermediate 
ranges  is  also  low  when  the  source  beam  is  diverted  by 
bottom  interaction  on  the  forward  trip.  However,  bottom 
interaction  on  the  return  trip,  which  must  occur  by  reci¬ 
procity,  often  causes  a  substantial  increase  in  overall  TL 
for  such  sites.  If  part  of  the  source  beam  is  rescattered  into 
CZ  propagation  from  these  intermediate  sites,  its  trajectory 
will  not  intersect  the  receiving  array.  Such  redistribution  of 


100  200  km 

FIG.  10.  Two-way  transmission  loss  map  of  TLf+TLb  for  observation 
R7E15.  A  general  correspondence  can  be  found  between  TLyH-TL*  min¬ 
ima  in  this  figure  and  backscatter  maxima  in  Fig.  7. 
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FIG.  11.  A  histogram  of  wide  area  TLy+TL*  for  observation  R7E15 
with  mean  level  191  dB  re:  1  m,  median  191  dB  re:  1  m,  and  standard 
deviation  14  dB.  Two-way  TL  values  corresponding  to  waterborne  sites 
comprise  the  lower  tail  of  this  distribution. 

acoustic  energy  is  only  of  significance  in  bistatic  observa¬ 
tion  geometries.  Due  to  uncertainty  in  geoacoustic  param¬ 
eters  of  the  seafloor,  the  level  of  modeled  TL  at  interme¬ 
diate  sites  insonified  after  previous  source-beam  bottom 
interaction  is  highly  approximate.  Since  the  angle  of  bot¬ 
tom  reflection  is  primarily  morphology  related,  the  loca¬ 
tion  of  these  intermediate  sites  is  more  dependable.  For 
sites  insonified  by  waterborne  paths,  which  are  more  cen¬ 
tral  to  the  analysis  of  this  paper,  knowledge  of  the  geoa¬ 
coustic  parameters  of  the  bottom  is  unnecessary  for  de¬ 
pendable  TL  modeling.  While  the  corresponding 
reverberation  map  in  Fig.  7  has  ambiguities  which  smear 
and  mirror  the  data  about  the  array’s  axis,  many  features 
are  closely  matched  by  the  two-way  TL  of  Fig.  10. 

The  associated  histogram  in  Fig.  11  shows  a  broad 
distribution  of  values  with  standard  deviation  significantly 
larger  than  that  of  the  backscatter  levels.  Sites  insonified  by 
waterborne  paths  correspond  to  the  lower  tail  of  the  dis¬ 
tribution.  Corresponding  high-level  backscatter  occupies 
the  upper  tail  of  the  distribution  shown  in  Fig.  8.  Like  the 
histogram  of  backscatter  level,  the  two-way  TL  distribu¬ 
tion  is  also  log  normal,  as  is  confirmed  by  the  0.98  normal¬ 
ized  correlation  coefficient  obtained  after  cross  correlation 
with  a  Gaussian  distribution  of  the  same  mean  and  vari¬ 
ance.  Again,  this  is  suggestive  of  an  intermittent  process, 
which  for  the  TL  case,  can  be  readily  attributed  to  the 
infrequent  intersection  of  the  source  beam  with  the  ocean 
bottom.  The  measurement  of  a  log-normal  distribution  for 
backscatter  level,  Fig.  8,  is  consistent  with  the  notion  that 
two-way  TL  drives  the  level  of  backscatter. 

To  facilitate  our  comparison  with  measured  backscat¬ 
ter,  we  incorporate  the  area  term  of  the  sonar  equation  by 
a  running  spatial  integration  of  the  antilog  of  negative  two- 
way  TL  as  windowed  by  the  spatially  varying  resolution 
rkr[3(6)  set  by  the  receiving  array’s  beamwidth  and  wave¬ 
form  bandwidth.  [We  substitute  the  antilog  of  TLy-j-TL* 
for  Ts(r>t )  in  Eq.  (4)  of  Ref.  10.]  This  incorporates  the 
area  term  of  the  sonar  equation  by  effectively  smearing  the 
TLf+TL*  map  to  the  resolution  of  the  backscatter  data 
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FIG.  12.  A  map  of  TLy+TL6— 10  log.4  for  observation  R7E15.  The 
area  term  is  incorporated  by  convolving  the  antilog  of  —  (TL^+TL&) 
with  the  spatially  varying  resolution  window  set  by  the  receiving  array’s 
beamwidth  and  waveform  bandwidth.  Right-left  ambiguity  is  not  in¬ 
cluded.  A  general  correspondence  can  be  found  between  minima  in  this 
figure  and  backscatter  maxima  in  Fig.  7. 

without  including  right-left  ambiguity,  as  is  shown  in  Fig. 
12.  The  new  map  is  referred  to  as  a  TLy+TLA— 10  log ^4 
map. 

It  is  now  evident  that  TL  selects  the  prominent  fea¬ 


tures  found  in  the  reverberation  map.  For  example,  at  long 
range  the  lineated  backscatter  maxima  at  site  a  is  precisely 
selected  by  an  extended  TL  minima,  confirming  our  previ¬ 
ous  deduction.  The  discrete  returns  associated  with  the 
waterborne  path  to  site  /J,  are  also  precisely  selected  by  TL. 
Also,  there  is  generally  higher  TLy+TLj,— 10  log  A  to  the 
west  and  east  of  the  observation  where  correspondingly 
low  backscatter  is  measured,  partly  due  to  the  north-south 
orientation  of  the  endfire  beams. 

While  a  broad  azimuthal  distribution  of  sites  in  the 
direct  path  have  distinct  TL  minima,  these  are  more  diffi¬ 
cult  to  compare  with  backscatter  due  to  the  excessive  over¬ 
lap  of  ambiguous  returns.  For  comparative  purposes,  we 
have  artificially  applied  right-left  ambiguity  to 
TL^+TL*— 10  log  ^4.  [We  substitute  the  antilog  of  the 
two-way  TL  for  f(r9t)  and  set  s(r,t)  =  1  in  Eq.  (5)  of  Ref. 
10.]  The  result  is  shown  in  Fig.  13.  A  great  similarity  can 
now  be  seen  between  backscatter  maxima  and  ambiguous 
TL^r+TL^— 10  log  A  minima  in  the  direct  path,  as  well  as 
at  longer  ranges  where  there  is  generally  a  pronounced 
TL^+TL^—  10  log  ^4  minima  for  every  prominent  back¬ 
scatter  maxima. 

To  quantify  this  comparison,  we  compute  the  normal¬ 
ized  cross  correlation  between  measured  backscatter  and 
modeled  TLyH-TL^— 10  log  ^4,  shown  in  Figs.  7  and  13, 
respectively.  We  convert  to  linear  units  which  are  propor¬ 
tional  to  root-mean-square  pressure  and  cross  correlate 
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FIG.  13.  Same  as  Fig.  13  except  right-left  ambiguity  is  now  included.  A  general  correspondence  can  be  found  between  minima  in  this  figure  and 
backscatter  maxima  in  Fig.  7.  See  Sec.  Ill  for  quantitative  cross  correlation  with  backscatter. 
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FIG.  14.  Reverberation  and  TL f  maps  for  observation  R7E10.  We  approximate  TL y+TL*  by  2TL y  for  R7E10,  and  again  incorporate  the  area  term  by 
a  running  integration.  Maps  of  2TLy— 10  log  A  with  and  without  ambiguity  appear  below.  See  Sec.  Ill  for  quantitative  cross  correlation  with  backscatter. 


jqR/20  with  j0-(TL/+TL<,-101og^)/20  The  rationale  for  this 

conversion  is  to  first  investigate  correlation  within  the 
physical  units  of  the  problem.  The  resulting  correlation 
coefficient  is  0.59.  However,  there  is  also  a  strong  physical 
motivation  for  correlating  these  quantities  directly  in  deci¬ 
bels,  as  they  are  presented  in  the  figures.  Besides  being 
linearly  related  in  the  sonar  equation,  our  analysis  shows 
that  backscatter  level  R  and  two-way  transmission  loss 
TLyr+TL*  have  normal  distributions.  Since  the  cross¬ 
correlation  function  is  a  second  moment  characterization, 
its  interpretation  is  unambiguous  for  variables  which  have 
normal  distributions.  That  is,  given  a  statistically  signifi¬ 
cant  number  of  samples,  a  high  correlation  unequivocally 
shows  dependence  while  a  low  correlation  shows  indepen¬ 
dence  for  Gaussian  variables.  Therefore,  the  normalized 
correlation  coefficient  of  0.63  obtained  after  cross  correlat¬ 
ing  R  with  (TL f+ TLb— 10  log  A )  is  not  only  quantitative, 
but  unambiguous  in  demonstrating  that  spatial  variations 
in  charted  backscatter  are  significantly  dependent  upon 
two-way  transmission  loss.  Refinements  in  charting  the 
backscatter  as  well  as  modeling  the  geoacoutic  parameters 
of  the  bottom  for  TL  computations  will  only  improve  the 
correlation. 

At  endfire,  however,  relative  backscatter  levels  are 
lower  than  are  predicted  by  modeled  TL f+TLb— 10  log  A. 
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This  is  most  likely  due  to  the  uniform  spatial  integration 
used  in  the  modeling.  Since  the  cross-range  extent  of  end- 
fire  sites  generally  exceeds  the  scale  of  geomorphological 
variation,  discussed  in  the  next  section,  significant  scatter¬ 
ing  strength  variation  is  likely  within  the  spatial  resolution 
of  endfire  sites.  This  could  easily  account  for  the  noted 
difference. 

We  generate  a  TLy  map  for  observation  R7E10,  as 
shown  in  Fig.  14.  Again,  there  is  a  strong  correspondence 
with  mapped  reverberation,  also  shown  in  Fig.  14  for 
R7E10.  In  particular,  the  highest  level  backscatter  features 
at  \\  CZ  are  precisely  selected  by  TL  minima  southwest  of 
the  observation  position.  This  site  is  in  the  general  vicinity 
of  the  elongated  feature  at  site  a  in  observation  R7E15. 
However,  the  magnitude  and  shape  of  these  returns  varies 
significantly  from  R7E15.  Since  both  observations  were 
within  15  km  CZ)  of  each  other,  we  conclude  that  small 
variations  in  measurement  position  can  lead  to  large  vari¬ 
ations  in  the  level  of  backscatter  at  a  given  site.  To  facili¬ 
tate  qualitative  comparison  with  backscatter,  we  approxi¬ 
mate  TLf+TLb  by  2TLy  for  R7E10,  and  again 
incorporate  the  area  term  by  a  running  integration.  The 
result  is  also  shown  in  Fig.  14  for  unambiguous  and  am¬ 
biguous  cases  respectively,  where  2TLy— 10  log  A  is  higher 
than  TL^+TL^— 101og,4  within  the  direct  path  due  to 
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FIG.  14.  Reverberation  and  TL f  maps  for  observation  R7E10.  We  approximate  TL^+TL*  by  2TL f  for  R7E10,  and  again  incorporate  the  area  term  by 
a  running  integration.  Maps  of  2TL f— 10  log  A  with  and  without  ambiguity  appear  below.  See  Sec.  Ill  for  quantitative  cross  correlation  with  backscatter. 


10*/2°  with  10-(TL/+TLf-101o^)/2°.  The  rationale  for  this 
conversion  is  to  first  investigate  correlation  within  the 
physical  units  of  the  problem.  The  resulting  correlation 
coefficient  is  0.59.  However,  there  is  also  a  strong  physical 
motivation  for  correlating  these  quantities  directly  in  deci¬ 
bels,  as  they  are  presented  in  the  figures.  Besides  being 
linearly  related  in  the  sonar  equation,  our  analysis  shows 
that  backscatter  level  R  and  two-way  transmission  loss 
TLy+TLj,  have  normal  distributions.  Since  the  cross¬ 
correlation  function  is  a  second  moment  characterization, 
its  interpretation  is  unambiguous  for  variables  which  have 
normal  distributions.  That  is,  given  a  statistically  signifi¬ 
cant  number  of  samples,  a  high  correlation  unequivocally 
shows  dependence  while  a  low  correlation  shows  indepen¬ 
dence  for  Gaussian  variables.  Therefore,  the  normalized 
correlation  coefficient  of  0.63  obtained  after  cross  correlat¬ 
ing  R  with  (TL y+TL*— 10  log  A )  is  not  only  quantitative, 
but  unambiguous  in  demonstrating  that  spatial  variations 
in  charted  backscatter  are  significantly  dependent  upon 
two-way  transmission  loss.  Refinements  in  charting  the 
backscatter  as  well  as  modeling  the  geoacoutic  parameters 
of  the  bottom  for  TL  computations  will  only  improve  the 
correlation. 

At  endfire,  however,  relative  backscatter  levels  are 
lower  than  are  predicted  by  modeled  TL y+TL^— 10  log  A. 
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This  is  most  likely  due  to  the  uniform  spatial  integration 
used  in  the  modeling.  Since  the  cross-range  extent  of  end- 
fire  sites  generally  exceeds  the  scale  of  geomorphological 
variation,  discussed  in  the  next  section,  significant  scatter¬ 
ing  strength  variation  is  likely  within  the  spatial  resolution 
of  endfire  sites.  This  could  easily  account  for  the  noted 
difference. 

We  generate  a  TLy  map  for  observation  R7E10,  as 
shown  in  Fig.  14.  Again,  there  is  a  strong  correspondence 
with  mapped  reverberation,  also  shown  in  Fig.  14  for 
R7E10.  In  particular,  the  highest  level  backscatter  features 
at  \\  CZ  are  precisely  selected  by  TL  minima  southwest  of 
the  observation  position.  This  site  is  in  the  general  vicinity 
of  the  elongated  feature  at  site  a  in  observation  R7E15. 
However,  the  magnitude  and  shape  of  these  returns  varies 
significantly  from  R7E15.  Since  both  observations  were 
within  15  km  CZ)  of  each  other,  we  conclude  that  small 
variations  in  measurement  position  can  lead  to  large  vari¬ 
ations  in  the  level  of  backscatter  at  a  given  site.  To  facili¬ 
tate  qualitative  comparison  with  backscatter,  we  approxi¬ 
mate  TLy+TL6  by  2TLy  for  R7E10,  and  again 
incorporate  the  area  term  by  a  running  integration.  The 
result  is  also  shown  in  Fig.  14  for  unambiguous  and  am¬ 
biguous  cases  respectively,  where  2TL y— 10  log  A  is  higher 
than  TLy+TL6— 10  log  ^4  within  the  direct  path  due  to 
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FIG.  15.  Lambertian  shaded  bathymetric  relief  for  a  light  source  centered 
above  the  region  at  infinity.  Site  a  corresponds  to  a  typical  abyssal  hill, 
and  /?  to  a  shallow  ledge  on  the  Mid-Atlantic  Ridge. 

the  above  approximation.  A  normalized  correlation  coeffi¬ 
cient  of  0.54  is  measured  between  backscatter  level  R  and 
the  combined  transmission  loss-area  terms  of  the  sonar 
equation  —  (TL y-f  TL^— 10  log  A )  for  observation 
R7E10.  Converting  to  linear  units  as  before,  the  correla¬ 
tion  coefficient  is  0.47.  These  results  are  consistent  with 
those  obtained  for  observation  R7E15. 

In  summary,  we  find  that  the  level  of  backscatter  is 
strongly  dependent  upon  TL.  Two-way  TL  selects  the  sites 
from  which  high-level  backscatter  is  returned,  with  back¬ 
scatter  maxima  corresponding  to  TL  minima.  Small  vari¬ 
ations  in  measurement  position,  less  than  |  CZ,  can  lead  to 
large  variations  in  the  level  of  backscatter  returned  from  a 
given  site  due  to  bathymetry  related  variations  in  TL. 
Larger  variations  in  measurement  position  can  also  lead  to 
spatial  variations  in  backscatter.  However,  these  may  sim¬ 
ply  be  due  to  shadow  zones  associated  with  CZ  propaga¬ 
tion,  and  are  not  necessarily  bathymetry  related. 

IV.  WESTERN  MID-ATLANTIC  RIDGE 
GEOMORPHOLOGY 

A.  Bathymetry  and  geomorphology 

Two  well-defined  corridors  of  excess  depth  run  across 
the  survey  area,  as  is  evident  in  Figs.  1  and  5.  These  are 
often  referred  to  as  “segment  valleys”  or  “fracture  zones” 
since  they  are  caused  by  separation  of  the  local  plate  seg¬ 
ment.  They  are  best  seen  in  Fig.  15  where  depth  perspec¬ 
tive  is  simulated  by  artificially  shading  the  bathymetry  as  a 
Lambertian  surface  with  a  light  source  at  infinity  directly 
above  center.  Darker  shading  indicates  steeper  slope.  The 
lower  segment  valley  runs  east-west  along  the  very  bottom 
of  the  image.  The  central  segment  valley  follows  a  more 
circuitous  southeast  to  northwest  route.  Ridge  boundaries 
found  on  the  northern  edge  of  the  central  segment  valley 


FIG.  16.  A  histogram  of  bathymetry  intersecting  the  circular  region  sur¬ 
veyed  by  backscatter  data  in  Fig.  7,  with  mean  depth  3785  m,  median 
depth  3833  m,  and  standard  deviation  467  m. 

are  believed  to  be  geologically  distinct  from  those  found  on 
southern  edge.  Northern  ridge  boundaries  are  known  as 
“inside  comers”  and  have  irregular  morphology.  Southern 
ridge  boundaries  are  known  as  “outside  comers”  and  have 
more  regularly  lineated  tendencies  typical  of  abyssal  hills. 
“Compared  to  outside  comers,  inside  comers  have  excess 
elevation,  high-residual  gravity  anomaly,  steep  and  often 
irregularly  oriented  slopes  created  by  multiple  high  angle 
faults...18.”19  Away  from  fracture  zones,  ridges  tend  to 
have  axes  parallel  to  the  MAR,  extending  from  North  to 
South.  They  also  tend  to  be  distributed  regularly  with  a 
spatial  period  of  about  5  km.  (The  faint  diagonal  striations 
throughout  the  shaded  bathymetry  at  about  20°  inclination 
from  horizontal  are  an  artifact  of  the  data  acquisition. 
Their  pattern  is  related  to  the  swaths  covered  by  the  re¬ 
search  vessel.)  A  histogram  of  the  depth  distribution  for 
the  entire  survey  area  is  presented  in  Fig.  16.  The  mean 
depth  is  within  a  standard  deviation  of  the  conjugate 
depth,  insuring  the  extensive  bottom  interaction  discussed 
in  previous  sections. 

B.  Directional  derivatives  of  bathymetry 

We  are  also  interested  in  the  relationship  between  re¬ 
verberation  and  seafloor  slope.  For  example,  does  high- 
level  backscatter  return  predominantly  from  steep  slopes? 
If  so,  is  it  biased  toward  slopes  facing  the  observation  po¬ 
sition?  To  study  this  we  have  computed  a  directional  de¬ 
rivative  (DD)  of  the  bathymetry  by  taking  the  dot  product 
of  the  gradient  vector  with  a  unit  vector  pointing  in  the 
direction  of  observation  R7E15  from  the  given  field  point. 
This  is  useful  since  it  gives  an  overview  of  the  vertical 
slopes  encountered  on  a  path  radially  outward  from  the 
source,  as  sound  travels.  A  spatial  map  of  DD  can  be  used 
to  detect  the  horizontal  orientation  of  ridges  with  vertical 
slopes  facing  the  observation  position.  The  result  is  shown 
in  Fig.  17.  Since  the  gradient  is  computed  by  finite  differ¬ 
ences,  the  DD  resolution  is  no  finer  than  the  200-m  grid 
interval  of  the  bathymetry. 

Here  we  identify  ridges  as  lineations  of  high  DD  mag¬ 
nitude.  For  a  given  ridge,  surfaces  facing  the  observation 
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FIG.  17.  A  directional  derivative  (DD)  of  the  bathymetry  equal  to  the 
dot  product  of  the  gradient  vector  at  each  point  in  the  field  with  the  unit 
vector  pointing  from  the  location  of  observation  R7E15.  This  gives  the 
slope  orientation  along  the  path  sound  waves  travel  from  the  source. 
Positive  values  indicate  slopes  facing  the  observation  (back  facing),  and 
negative  values  indicate  slopes  facing  away  (forward  facing).  Slopes  are 
given  by  tan  <p,  where  <p  is  the  slope  angle  from  horizontal  along  the  DD. 

point  have  positive  DD,  surfaces  facing  away  have  negative 
DD,  and  surfaces  normal  have  zero  DD.  Since  the  ridge 
axes  generally  run  north-south,  we  often  find  north-south 
lineations  of  high  DD  magnitude  east  and  west  of  the  ob¬ 
servation  point  corresponding  to  abyssal  hills.  To  the  north 
and  south,  ridge  axes  are  generally  parallel  to  the  observa¬ 
tion  direction  where  we  find  DD  maxima  diffusely  distrib¬ 
uted  in  speckled  patterns  corresponding  to  ridge  edges  and 
crests.  These  patterns  are  better  defined  at  segment  valley 
walls  along  the  fracture  zone  central  to  the  observation, 
and  are  most  coherent  for  the  irregularly  shaped  inside 
comers  which  give  much  higher  DD  values  than  outside 
comers,  as  is  anticipated.  Abyssal  and  segment  valley 
floors  tend  toward  zero  slope. 

A  histogram  of  the  DD  distribution  for  the  survey  area 
is  presented  in  Fig.  18  where  a  sharply  peaked  and  sym¬ 
metric  distribution  with  8°  variance  about  zero  slope  is 
found.  Slopes  facing  the  observation  position  are  as  likely 
to  be  found  along  a  propagation  radial  as  slopes  facing 
away.  However,  tails  of  the  distribution  beyond  ±25°  are 
not  accurate.  This  is  due  to  limitations  of  the  measurement 
system  and  coarseness  of  the  sampling  interval. 

To  obtain  a  quantitative  estimate  of  the  scales  over 
which  bathymetry  varies,  we  have  computed  the  wave- 
number  spectrum,  shown  in  Fig.  19.  The  dominant  spec¬ 
tral  component  is  in  the  low  wave-number  regime  corre¬ 
sponding  to  length  scales  greater  than  about  5  km. 
Variations  within  the  array’s  range  resolution  ( 1.5  km)  are 
less  significant  on  average.  However,  from  the  DD  map  in 
Fig.  17,  steep  slopes  do  exist  within  these  smaller  scales. 
While  they  are  less  frequently  encountered,  steep  slopes 
associated  with  lineated  ridges  are  more  likely  to  terminate 
shallow-angle  propagation,  typical  at  long  ranges,  than 
level  areas  which  tend  to  occur  in  valleys  below  the  con¬ 
jugate  depth  contour.  High  wave  number  spikes  of  roughly 


FIG.  18.  A  histogram  of  DD  intersecting  the  circular  region  surveyed  by 
backscatter  data  in  Fig.  7,  with  mean  slope  0.20°,  median  0.23°  and  stan¬ 
dard  deviation  8.0°  from  horizontal.  Slopes  near  zero  predominate,  and 
neither  slopes  facing  or  opposing  the  observation  are  favored. 

20  dB  occur  for  spatial  frequencies  greater  than  1  km”1. 
These  are  an  artifact  of  registration  error  along  neighbor¬ 
ing  measurement  swaths  and  do  not  reliably  correspond  to 
actual  bathymetric  variations. 

V.  COMPARISONS  WITH  BATHYMETRY 

To  correlate  backscatter  with  geomorphology,  a  color 
reverberation  map  is  projected  onto  the  black  and  white 
shaded  bathymetry  relief.  The  resulting  image  is  presented 
in  Fig.  20  for  observation  R7E15.  To  insure  that  bathy¬ 
metric  variations  are  not  obscured  in  the  relief,  only  high 
level  backscatter  is  presented. 

While  cross-range  smearing  and  right-left  ambiguity 
make  the  image  somewhat  confusing,  close  inspection  re- 


FIG.  19.  A  two-dimensional  wave  number  spectrum  of  hydrosweep 
bathymetry  over  a  102.4x102.4  km2  spatial  window  within  the  region 
shown  in  Fig.  1.  The  dominant  spectral  components  occur  for  length 
scales  greater  than  about  5  km.  High  wave-number  spikes  of  roughly  20 
dB  occurring  for  spatial  frequencies  greater  than  1  km-1  do  not  corre¬ 
spond  to  actual  bathymetric  variations.  These  are  an  artifact  of  registra¬ 
tion  error  along  neighboring  measurement  swaths  of  the  bathymetry. 
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FIG.  20.  Color  backscatter  levels  from  observation  R7E15  projected  onto  corresponding  Lambertian  shaded  bathymetry,  see  Fig.  16.  Right-left 
ambiguity  has  not  been  eliminated.  Higher  level  backscatter  returns  often  resemble  the  bathymetric  features  to  which  they  are  mapped.  Only  prominent 
backscatter  is  shown  to  avoid  obscuring  bathymetric  variations. 


veals  that  many  high-level  backscatter  features  closely 
match  corresponding  bathymetry.  However,  we  do  not 
base  our  distinction  between  true  and  virtual  scattering 
sites  upon  visual  correlation.  While  this  method  may  often 
be  accurate,  as  well  as  intuitive,  the  probability  of  false 
correlation  is  also  high  given  cross-range  resolution  which 
is  often  only  slightly  less  than  the  separation  between 
ridges. 

Instead,  we  resolve  right-left  ambiguity  by  comparing 
TLy+TL*— 10  log  ^4  at  ambiguous  scattering  sites.  If  the 
difference  in  TL y+TL^— 10  log  A  is  greater  than  a  chosen 
threshold,  we  associate  the  true  scattering  site  with  the 
place  of  lower  TLy+TL^— 10  log  ^4  and  map  the  corre¬ 
sponding  backscatter  levels  to  this  site.  If  the  difference  is 
less  than  this  threshold  we  do  not  make  a  distinction,  and 
project  backscatter  level  onto  neither  site.  We  choose  a 
threshold  of  10  dB  which  constitutes  a  significant  enough 
difference  in  TLy+TL*— 10  log  A  to  rule  out  false  identi¬ 
fication  in  most  cases. 

The  resulting  unambiguous  backscatter  levels  are  pro¬ 
jected  onto  shaded  bathymetry  and  presented  in  Fig.  21  for 
observation  R7E15.  Again,  only  high-level  backscatter 
( >  70  dB  re:  1  /xPa)  is  presented  to  insure  that  bathymetric 
variations  are  not  obscured.  We  first  note  that  folding  the 
figure  about  the  array’s  axis  yields  the  same  prominent 
returns  depicted  in  the  ambiguous  Figs.  7  and  20.  In  other 
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words,  little  information  about  high-level  backscatter  is 
lost  by  the  choice  of  a  10-dB  threshold.  On  the  other  hand, 
a  general  correspondence  between  prominent  bathymetry 
and  high-level  backscatter  is  now  evident.  In  particular, 
the  backscatter  lineation  at  site  a  discussed  in  previous 
sections  corresponds  to  an  abyssal  hill  to  the  southwest  of 
the  observation.  The  well-defined  pattern  of  high-level 
backscatter  found  in  the  direct  path  traces  the  ridge  con¬ 
tours  surrounding  the  segment  valley  which  runs  central  to 
the  observation  position. 

The  relationship  between  prominent  backscatter  and 
the  ridge  structure  is  further  elucidated  in  Fig.  22  by  over¬ 
lay  on  a  DD  map.  Here  prominent  returns  can  be  more 
easily  associated  with  specific  ridges.  For  example,  this  is 
the  case  within  the  direct  path,  as  well  as  at  long-range 
sites  a  and  /?.  However,  it  is  also  evident  that  the  resolution 
window  is  frequently  wider  than  scarps  within  the  ridges. 
That  is,  returns  are  typically  mapped  to  areas  larger  than 
the  width  of  associated  scarps.  As  a  result,  quantitative 
comparisons  between  slope  orientation  and  backscatter  are 
not  possible  with  the  given  measurements,  and  are  left  for 
future  analysis  with  higher  resolution  data.  (Preliminary 
results  from  high  resolution  measurements  made  during 
the  main  acoustics  experiment  show  a  strong  relationship 
between  backfacing  scarps  and  backscatter.20) 

However,  an  association  between  backfacing  scarps 
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FIG.  21.  Same  as  Fig.  20  except  ambiguity  has  been  eliminated  with  a  TL-based  environmental  symmetry  breaking  technique.  Only  prominent 
backscatter  is  shown  to  avoid  obscuring  bathymetric  variations.  Highest  level  backscatter  is  found  where  waterborne  paths  are  available.  Well-defined 
patterns  in  backscatter  often  correspond  to  bottom  morphology,  following  ridges  protruding  above  the  conjugate  depth  contour. 


FIG.  22.  White  contours  of  unambiguous  backscatter  exceeding  80  dB  re: 
1  /iPa  are  projected  onto  the  color  map  of  DD  shown  in  Fig.  17.  Well- 
defined  patterns  in  backscatter  often  follow  ridges  with  extended  back- 
facing  scarps.  Backscatter  resolution  for  the  given  signal  is  generally  too 
poor  resolve  the  scarps. 


and  prominent  backscatter  can  still  be  made  in  the  present 
study  with  the  aid  of  the  TL  modeling  presented  in  Secs.  I 
and  III.  Specifically,  the  resemblance  between  high-level 
backscatter  and  bathymetry  can  be  largely  attributed  to 
spatial  variations  in  TL.  We  focus  on  the  numerous  sites 
which  are  insonified  by  waterborne  paths  and  therefore 
return  the  most  prominent  backscatter.  This  is  because 
there  is  little  question  regarding  the  wide-angle  PEs  accu¬ 
racy  in  modeling  TL  along  waterborne  paths  for  a  range- 
independent  water  column  in  the  deep  ocean.  For  example, 
compare  the  TL  maps  in  Figs.  10  and  12  with  the  bathy¬ 
metric  maps  shown  in  Figs.  1,  5,  15,  and  17.  Within  the 
direct-path,  segment  valley  ridges  facing  observation 
R7E15  are  broadly  outlined  by  TL  minima  along  the  con¬ 
jugate  depth  contour.  And,  lineations  of  minimum  TL  oc¬ 
cur  at  1^  CZ  ranges  on  backfacing  scarps  for  ridges  which 
protrude  above  the  conjugate  depth  contour  over  a  wide 
azimuth,  as  is  the  case  at  site  a .  In  general,  ridges  which 
cross  the  propagation  path  are  more  likely  to  show  up  as 
TL  lineations  since  they  span  broader  azimuthal  sectors 
than  those  running  parallel  to  the  propagation  path.  As  a 
consequence  of  the  predominantly  north-south  running 
ridge  axes  of  the  western  MAR,  ridges  east  and  west  of  the 
observation  are  more  likely  to  yield  backscatter  lineations. 
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FIG.  23.  Scattering  strength  over  a  wide-area  for  observation  R7E15.  The 
estimate  is  based  upon  the  unambiguous  backscatter  levels,  see  Fig.  21 
and  caption,  and  the  unambiguous  TLy+TL*— 10  log  A  levels  presented 
in  Fig.  12  for  a  233  dB  re:  1  /iPa  and  1-m  source  strength. 

VI.  SCATTERING  STRENGTH 

We  estimate  scattering  strength  5  over  a  wide  area  by 
adding  unambiguous  backscatter  level  R  to  modeled 
TLy+TL^— 10  log  ^4  and  subtracting  the  source  strength 
JV=233  dB  re:  1  /xPa  and  1  m.  The  result  appears  in  Fig. 
23  for  observation  R7E15  where  we  find  uniformly  low 
scattering  strength  values.  (Here  we  have  used  the  full 
range  of  unambiguous  backscatter  levels  as  opposed  to  just 
the  prominent  values  overlain  on  relief  bathymetry  in  Fig. 
21.)  A  histogram  is  shown  in  Fig.  24  where  scattering 
strength  is  distributed  narrowly  about  the  median,  roughly 
following  the  anticipated  log-normal  distribution  found  in 
backscatter  level  and  transmission  loss.  A  0.96  normalized 
correlation  coefficient  is  obtained  after  cross  correlation 
with  a  Gaussian  distribution  of  the  same  mean  and  vari¬ 
ance.  The  low  variance  and  generally  uniform  spatial  dis- 


FIG.  24.  A  histogram  of  the  scattering  strength  values  at  sites  presented 
in  Fig.  23,  with  mean  —47  dB,  median  —47  dB  and  standard  deviation  7 
dB.  Since  the  signal  used  cannot  resolve  distinct  bathymetric  features 
which  may  dominate  the  backscatter  process,  these  statistics  are  primarily 
useful  as  a  lower  bound  for  the  survey  area. 
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tribution  of  scattering  strength  is  consistent  with  the  low 
resolution  of  the  measurement.  That  is,  contributions  from 
elemental  scatterers  are  well  mixed  in  the  long-range,  low- 
resolution  backscatter  returns. 

The  scattering  strength  values  presented  are  useful  as  a 
lower  bound  for  the  region.  This  is  because  the  backscatter 
data  presented  cannot  resolve  the  dominant  contributors, 
as  is  suggested  in  the  previous  section.  Since  these  contrib¬ 
utors  are  under-resolved,  their  scattering  strength  is  likely 
underestimated.  (As  a  point  of  interest,  we  note  that  the 
lower  bound  measured,  47  ±4  dB  is  roughly  consistent 
with  Lambert’s  law  5=20  log(sin  i/0  ~(x  and  the  Macken¬ 
zie  coefficient  fi=  —  27  dB  (Ref.  21 )  for  the  predominantly 
shallow  grazing  angles,  tp<  10,  of  this  survey.  However,  we 
believe  that  this  result  is  merely  coincidental.  Recent  work 
with  high-resolution  FM  signals  in  the  same  area  indicate 
that  actual  scattering  strengths  are  significantly 
higher.16-22) 

VII.  CONCLUSIONS 

Two  primary  conclusions  are  drawn  from  this  initial 
analysis  of  long-range  backscatter  in  the  western  MAR: 
(a)  there  is  a  high  correlation  between  prominent  back¬ 
scatter  and  detailed  bottom  morphology;  (b)  the  range  and 
azimuth  of  prominent  backscatter  can  be  determined  by 
propagation  modeling  if  high-resolution  bathymetry  is 
available.  In  general,  measured  backscatter  maxima  corre¬ 
spond  to  modeled  transmission  loss  (TL)  minima  when 
charted  to  respective  locations  on  the  sea  floor.  These  spa¬ 
tial  extrema  are  most  pronounced  in  the  direct  path  and  at 
n+\  CZ  ranges  when  propagation  is  not  blocked  at  shorter 
range  for  the  given  azimuth,  where  n= 0,1  for  the  present 
study.  At  these  ranges,  backscatter  maxima  typically  fol¬ 
low  extended  ridges  along  the  local  conjugate  depth  con¬ 
tour.  This  is  partially  explained  by  propagation  modeling 
which  shows  that  forward  propagation  paths  tend  to  in- 
sonify  the  backfacing  scarps  along  ridges  which  protrude 
above  the  conjugate  depth  contour.  This  leads  to  the  cor¬ 
responding  TL  minima  which  apparently  play  a  fundamen¬ 
tal  role  in  driving  the  level  of  backscatter.  Whether  the 
prominence  of  these  returns  is  also  related  to  scattering 
properties  of  the  associated  ridges  and  their  orientation  is 
still  uncertain.  This  uncertainty  arises  because  the  resolu¬ 
tion  of  the  measurements  presented  in  this  study  is  often 
too  poor  to  adequately  resolve  the  scarps,  let  alone  more 
elementary  contributors  to  the  scattering  process  which 
may  reside  within  the  scarps. 

Although  we  do  not  yet  fully  understand  the  scattering 
process,  we  have  learned  a  great  deal  about  reverberation. 
Our  analysis  has  provided  a  simple  method  of  forecasting 
the  range  and  azimuth  of  prominent  backscatter  in  the 
MAR.  Either  by  a  crude  tracing  of  CZ  circles  on  an  excess 
depth  map,  or  a  more  accurate  modeling  of  transmission 
loss  with  a  full-wave  propagation  model  such  as  the  wide- 
angle  PE.  Both  of  the  above  methods  were  used  in  design¬ 
ing  the  main  acoustics  experiment  of  the  SRP.  The  prelim¬ 
inary  results  of  this  more  recent  experiment  show  that  both 
methods  are  extremely  dependable  in  forecasting  promi¬ 
nent  backscatter  in  the  MAR.20,23 
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We  have  also  estimated  scattering  strength  over  a  wide 
area  and  have  obtained  a  lower  bound  of  —  47  ±4  dB  for 
the  survey.  Our  estimates  provide  a  lower  bound  because 
well-defined  morphological  features,  which  apparently 
play  a  major  role  in  the  backscatter  process,  are  under- 
resolved  by  the  cw  measurements  used  in  this  study.  More 
recent  scattering  strength  estimates  made  with  much 
higher  resolution  FM  data  have  yielded  significantly  higher 
values  for  the  same  region.22  Nevertheless,  the  estimates 
presented  here  show  a  low  variance  and  generally  uniform 
spatial  distribution.  This  is  consistent  with  the  high  cross 
correlation  measured  between  long-range  backscatter  and 
TL  in  this  study.  Evidently,  distinct  scattering  contributors 
are  so  well  mixed  in  the  large  areas  integrated  by  long- 
range  returns  that  the  most  significant  variations  in  back¬ 
scatter  are  generally  due  to  propagation.  This  further  sup¬ 
ports  the  idea  of  backscatter  forecasting  by  propagation 
modeling.  It  also  suggests  that  the  possibility  of  remote 
bathymetric  survey  by  backscatter  measurement5,6,24  is 
practical.  However,  a  more  negative  implication  is  that  any 
of  a  number  of  scattering  models  will  be  able  to  match 
long-range  MAR  backscatter  data,  regardless  of  whether 
they  are  correct,  so  long  as  a  sufficient  propagation  model 
is  used. 

Investigations  into  the  fundamental  scattering  process 
are  better  suited  to  the  site-specific  Main  Acoustics  Cruise 
of  July  1993.  This  bistatic  experiment  was  designed  specif¬ 
ically  around  three  sites  detailed  in  the  hydrosweep 
bathymetry.  High-resolution  monostatic  and  bistatic  rever¬ 
beration  measurements  were  made  at  various  ranges  and 
orientations  from  known  scarps.  Short-range  measure¬ 
ments  should  enable  resolution  of  distinct  small-scale  con¬ 
tributors  to  the  scattering  process.  While,  long-range  mea¬ 
surements  will  show  how  these  accrue  in  distant 
reverberation.  Vertical  arrays  were  also  deployed  before 
select  backfacing  scarps  to  explore  the  vertical  structure  of 
the  scattered  field  near  the  scatterers.  These  data  are  to  be 
compared  with  hydrosweep  as  well  as  finer  scale  bathym¬ 
etry  and  geological  data  recently  collected  to  further  elu¬ 
cidate  the  scattering  process.  In  addition,  tow-ship  tracks 
in  the  Main  Acoustics  Experiment  were  designed  with  am¬ 
biguity  removal  as  an  essential  theme.  Star  patterns  rather 
than  polygons  were  used  extensively  to  provide  redundant 
measurements  at  differing  array  headings  at  a  point  of  con¬ 
vergence.  This  was  necessary  because  the  present  study  has 
shown  that  small  differences  in  measurement  position,  as 
small  as  CZ/4,  can  lead  to  drastic  differences  in  backscat¬ 
ter  due  to  bathymetry  related  variations  in  TL.  Preliminary 
results  from  the  ’93  cruise  indicate  that  the  environmental 
symmetry  breaking  technique  for  ambiguity  removal  pre¬ 
sented  here  is  consistent  with  ambiguity  removal  by  com¬ 
parison  of  data  from  multiple  adjacent  array  headings.20  A 
detailed  report  on  these  results  is  in  preparation.23 
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